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ABSTRACT 

If star formation proceeds by thermal fragmentation and the subsequent gravitational collapse of the 
individual fragments, how is it possible to form fragments massive enough for O and B stars in a typical 
star-forming molecular cloud where the Jeans mass is about 1 M Q at the typical densities (10 4 cm -3 ) 
and temperatures (10 K)? We test the hypothesis that a first generation of low-mass stars may heat 
the gas enough that subsequent thermal fragmentation results in fragments >10M Q , sufficient to form 
B stars. We combine ATCA and SMA observations of the massive star-forming region G8. 68— 0.37 
with radiative transfer modeling to derive the present-day conditions in the region and use this to 
infer the conditions in the past, at the time of core formation. Assuming the current mass/separation 
of the observed cores equals the fragmentation Jeans mass/length and the region's average density 
has not changed, requires the gas temperature to have been 100 K at the time of fragmentation. The 
postulated first-generation of low-mass stars would still be around today, but the number required 
to heat the cloud exceeds the limits imposed by the observations. Several lines of evidence suggest 
the observed cores in the region should eventually form O stars yet none have sufficient raw material. 
Even if feedback may have suppressed fragmentation, it was not sufficient to halt it to this extent. 
To develop into O stars, the cores must obtain additional mass from outside their observationally 
defined boundaries. The observations suggest they are currently fed via infall from the very massive 
reservoir (^1500 Mq) of gas in the larger pc scale cloud around the star-forming cores. This suggests 
that massive stars do not form in the collapse of individual massive fragments, but rather in smaller 
fragments that themselves continue to gain mass by accretion from larger scales. 
Subject headings: stars: formation — submillimeter — stars: winds, outflows — ISM: clouds — stars: 
early-type — ISM: evolution 



1. INTRODUCTION 

Massive stars play an influential role in shaping the 
Universe. Observation s show that the su bstantial ma- 
jority form in clusters (jLada fc L ada 2003:) yet the phys- 
ical processes governing the fragmentation and collapse 
of their natal molecular cloud, a crucial step in determin- 
ing important parameters such as the number of massive 
stars and their final stellar mass, remains an unsolved 
problem. Infrared dark clouds, thought to be examples 
of the birth sites of massive clusters, are observed to 
have temperatures of 10 — 20 K and contain many hun- 
dreds to thousands of so l ar masses of gas dPillai et al. 
2006] iRagan et all 12001 IRathborne et al l 120061: iSwift 
2003: iZhang et alJ 1200ft IRathborne et all 120101 ) . It re- 
mains an open question whether these initial conditions 
can produce fragments large enough to form high-mass 
stars through direct collapse, without the cores them- 
selves sub- fragmenting first. Radiative feedback (heat- 
ing) from embedded, low mass protocluster members has 
been proposed as a mechanism to delay f ragmentation by 
chang ing the effective equation of state (jKrumhol z et al.l 
120071 ). In this paper we consider a more direct mecha- 
nism to suppress the sub-fragmentation. Through obser- 
vations of a very young high-mass protocluster, we inves- 
tigate whether heating from a first generation of low mass 
stars can raise the cloud temperature enough that subse- 
quent thermal fragmentation can produce fragments with 

1 Contact email: slongmore@cfa.harvard.edu 



sufficient mass to form high-mass stars. 

1.1. The IRAS 18032-2137 complex 

Figure Q] shows the IRAS 18032-2137 star forming 
complex, comprised of 3 distinct regions separated by 
a few arcminutes on the sky. The most evolved of 
these is the stellar cl uster, BDS2003-3 (a.T2nnn=18:06:15, 
<5j2ooo = -21:37:30, iBica et al.l l200l . The high vi- 
sual extinction, association with both near-IR nebulosity 
(|Longmore fc Burton! 120091) and radio continuum emis- 
sion fG8.662-0.342. iBecker et al.l 119941 ) make this most 
likely a massive, heavily embedded cluster. One ar- 
cminute east of the cluster lies the ultr a-comp act HII 
region G8.67 - 0.36 (jWood fc Churchwel]lll98^_ coinci- 
dent with a compact sub- mm cont inuum peak (|Hill et al.l 
120051 120061: iThompson et al.ll2006D . One arcminute (pro- 
jected separation of 1.4 pc) north-east along the same 
sub-mm continuum filament lies a second sub-mm con- 
tinuum peak, G8.68 - 0.37. Both G8.67 - 0.36 and 
G8.68 - 0.37 are associated with H 2 Q, class II CH3OH 
and OH maser emission (Hofner fc Churchwelll 1996; 



Wals h et al.l 19 98: Forstcr fc Caswel|ll98ftlCaswelllll998l: 
Val'tts et al]|2000ft . iHill et al.l (|2010l) show G8.68 - 0.37 
is massive (^1.5xl0 3 M Q ) and has a high luminosity 
(~1O 4 L ) yet deep cm-continuum observations reveal 
no com pact (Michele Pestallo zi priv. comm.) or ex- 
tended (|Longmore et al. | 120091) free-free e mission. At 
a distance of 4.8 kpc ( Purcell et al.l [20061 ) . the 3cr up- 
per limit of 1.7mJy at 3.6 cm (Michele Pestallozi priv. 
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comm.) corresponds to a Lyman-continuum photon rate 
of 1.65xl0 46 s _1 , implying no stars earlier than B0-B1 
are present in the cluster. All of this points to the na- 
ture of G8. 68 — 0.37 as a very young massive star forming 
region, prior to the formation of ultra-compact or hyper- 
compact HII regions. 

NH3 observations of G8. 68— 0.37 show: (i) an extended 
cold component seen in the NH3(1,1) & (2,2) transi- 
tions with a morphology similar to the sub-mm contin- 
uum emission; (ii) a warmer component in NH3(4,4) & 
(5,5), unresolved at 8" at the peak of the sub-mm con- 
tinuum and methanol maser emission (IPillai et alJl2007t 
iLongmore et alJl2007t ). This suggests G8.68 - 0.37 is in- 
ternally heated by young (proto)stars of at least several 
M©. Strong infall profiles are seen in 3mm Mopra spec- 
tra o f HCO +, HNC and e ven 1 3 CO (jPurcell et alj [20061 
120091) . while lHariu et al.1 ([19981 ) report strong SiO emis- 
sion indicative of shocks caused by outflows - further 
evidence that star formation is already underway in the 
protocluster. 

In summary, G8.68 — 0.37 is a relatively isolated mas- 
sive star forming core with a symmetric and centrally 
peaked dust profile. It is still in the earliest stages of 
forming a massive protocluster and has already begun 
to significantly heat the gas in the proto-cluster center. 
As such it appears a good candidate region for testing 
whether energetic feedback can act sufficiently quickly 
to suppress thermal fragmentation. 

2. OBSERVATIONS AND DATA REDUCTION 

The data were taken with the Submillimeter Arraj0 
(SMA) between 2007 September 2 nd and 2008 September 
17 th in three individual tracks between 217 and 279 GHz 
in the sub-compact, compact and extended array config- 
urations. At these frequencies, the SMA primary beam 
(field-of-view) is ~ 45" — 58", sufficient to cover the ex- 
tent of the single-dish continuum emission (see Figure [1} 
in a single pointing centered at aj2ooo =18:06:23.47, 
^72000 = —21:37:7.6. For each track, this sky position 
was observed for periods of 10-15 minutes on-source, in- 
terspersed with observations of two bright, nearby cal- 
ibrators (1733-130 k 1911-201). Two bandpass cali- 
brators, 3c279 and 3c454.3, were observed for approxi- 
mately an hour at the start and end of each observation. 
At least one primary flux calibrator (Uranus, Neptune, 
Titan) was observed for each track. The absolute flux 
scale is estimated to be accurate to ^15%. In all cases 
the weather was very stable and the resulting ampli- 
tude/phase stability was good. The data were calibrated 
using the MIR IDL packag43 and exported to Miriad to 
be imaged and cleaned. A zeroeth order polynomial was 
fitted to the line-free channels and subtracted from the 
visibilities to separate the line and continuum emission, 
which were imaged separately. Table [T] lists the observ- 
ing setup, continuum sensitivity and resolution for each 
of the tracks. The results are shown in § 13.11 

3. RESULTS 

2 The Submillimeter Array is a joint project between the Smith- 
sonian Astrophysical Observatory and the Academia Sinica Insti- 
tute of Astronomy and Astrophysics and is funded by the S mith- 
sonian Institution and the Academia Sinica (Ho et al. 2003). 

3 http:/ /www. cfa.harvard.edu/~cqi/mircook. html 



3.1. SMA Data 

In the sub-compact and compact configuration data, 
the continuum morphology is similar to that seen in the 
single-dish observations - a single component, centrally 
peaked towards the methanol maser position. However, 
as shown in Figure [2j the SMA extended configura- 
tion observations resolve the emission into 3 components, 
named MM1 to MM3 in order of peak intensity. These 
are separated by between 1.4"and 2.2", corresponding to 
projected distances of ~6200AU and ^9700 AU, respec- 
tively. Table [5] gives the properties of the mm continuum 
detections. 

Figures [3] & |4] show the 230 GHz lower and upper side- 
band spectra, respectively, towards the peak of the con- 
tinuum emission from the compact configuration data. In 
addition to the 13 CO (2-»l) and 12 C0(2->1) emission, 
the spectra show several more complex molecules (e.g. 
CH3OH & CH3CN) confirming the temperature is suffi- 
ciently high for these molecules to have evaporated off the 
dust grains and into the gas phase. However, the spec- 
tra do not display the rich inventory of molecular l ines 
seen towards hot molecular cores (see lCesaronll2005l for 
a review) suggesting this region is an intermediate evo- 
lutionary stage between the cold and hot core stages. A 
similar result is found for the extended configuration ob- 
servations which were tuned to 217 GHz (instead of 230 
GHz). Figure [5] shows the 12 CO (2— >1) channel maps in 
which a bipolar morphology is clearly seen at opposite 
sides of the systemic Vlsr (37.2 kins" 1 ), indicative of 
a molecular outflow. The properties of this outflow are 
discussed in more detail in § 14.21 

The emission from more complex species detected in 
the 230 GHz compact configuration is generally unre- 
solved within the 3"x2" beam and encompasses MM1 to 
MM3. This emission is also coincident with the NH3(4,4) 
& (5,5) and methanol maser emission. The spectral lines 
are typically well fit as a single Gaussian at the VlsRj 
with linewidth ^5kms _1 and no sign of any coherent, 
large-scale velocity structure. 

3.2. Spitzer Data 

Figure [21 shows a 3-color GLIMPSE image 
(jBeniamin et all I2003D at 3.6, 4.5 and 8.0 /iin to- 
wards the region taken with IRAC on the Spitzer Space 
Telescope. There is a clear bipolar nebulosity in the 
image, prominent in all four IRAC bands. The black 
contours show the positions of the three mm-continuum 
cores detected with the SMA relative to the IRAC 
emission. Of the three mm cores, the nebulosity is most 
closely associated with MM1. In the longer wavelengths 
(24 & 70 urn) and lower r esolution (7.5 - 22") MIPSGAL 
data dCarev et al.l 12009ft . only a single, bright (1 Jy 
integrated flux at 24 (im) emission peak is seen with 
MIPS. 

4. ANALYSIS 

4.1. Total Luminosity of G8. 68-0. 37 

While luminosity estimates for this region are available 
from the literature, access to recent Spitzer and deep 
near-IR data not available to previous authors enables 
us to provide tighter constraints. We used the Spitzer 
24 & 70 /xm, SCUBA 450 & 850 ^m and SIMBA 1.3mm 
(Hill et al. 2005) data to estimate the total luminos- 
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ity of the system. Deep, near-IR imaging show this re- 
gion is infrared-dark down to 18th magnitude at 2.2/im 
(|Longmore fc B urton 2009) , giving confidence that most 
of the luminosity is emitted at the longer wavelengths we 
used to construct the spectral energy distribution. Based 
on the similar size and morphology of the emission from 
24/im to 1.3mm, we used a circular aperture of radius 
15" to derive the total flux at each wavelength. We fit 
the resulting spectral energy distribution as a greybody 
and, integrating under the resultant fit, estimate the pro- 
tocluster luminosity to be ^1.9xl0 4 L Q . The SMA data 
were not used in the luminosity estimate as comparison 
with the SIMBA data (see § I4.3.1[) revealed a large frac- 
tion of the total flux has been spatially filtered by the 
interferometer. 

4.2. Outflow Properties 

As shown in Figure [5J both the 230 GHz compact- 
configuration, CO SMA data and Spitzer data reveal a 
bipolar outflow in an approximately NE-SW orientation. 
The higher angular resolution, extended-configuration, 
217GHz continuum SMA data shows core MM1 lies di- 
rectly along the projected outflow axis of both the CO 
and IRAC nebulosity making it most likely to be driv- 
ing the outflow. Figure [5] shows channel maps of the 
12 CO (2— >T) emission. Extended emission from the ambi- 
ent cloud is filtered out by the interferometer (our short- 
est baseline corresponds to a spatial filtering of emission 
more extended than 30") close to the LSR velocity of 
37.2kms- 1 (Purcell et al. 2006, 2009). To determine the 
outflow properties we first separated the outflow into its 
two component lobes using a 3 sigma cutoff in the 12 CO 
data cube at velocities 29.3-36.3 km/s (blue) and 45-74 
km/s (red). Velocities closer to the LSR velocity were 
avoided to prevent contamination from the ambient low 
density gas. 

The outflow mass was derived following IScoville et al.l 
(fl98ll . assuming a CO excitation temperature of 30 K. 
We defined a blue and red polygon covering the spatial 
extent of the outflow emission shown in Figure [3J and 
derived an average spectrum for both 12 CO & 13 CO. 
Where 13 CO was detected, we estimated the 12 CO op- 
tical depth by assuming the 13 CO emission to be opti- 
cally thin with an abundance ratio to 12 CO of 89 (see 
I Wilson fc Matteuccilfl992l) . The 12 CO optical depth was 
then used to compute the CO column density. For veloc- 
ities with no 13 CO emission, optically thin 12 CO emis- 
sion was assumed. We thus derive a total outflow mass 
of ~6M G . 

Based on the distance between the edge of the lobes 
through the projected center, the outflow has a dynami- 
cal time of ~1.3xl0 4 years. The derived outflow rate is 
therefore ~4.4xl0~ 4 M Q per year. We emphasize that 
due to missing flux (see Figure [5]) and optical depth ef- 
fects, the mass and outflow rate could be underestimated. 
Still, for the luminosity of the protocluster, the outflow 
properties are in good agreement with those seen to- 
war dsothCThjghj^ regions (see for exam- 
ple. IShepherd fc Churchwelll [1996 : IBeuther et al.l [200l 
iZhang et al.ll2001l 120051: lOiu et alll2009l ). 

Outflow studies in high mass star forming regions are 
often hampered by the complexity of the observed emis- 
sion. In most cases one sees multiple outflows, many 



with unknown driving sources, which may be difficult to 
disentangle. While there are potentially signs of addi- 
tional outflows in Figure [5] (e.g. the SE clump at offset 
5", -10" and 31.5 to 34.5kms" 1 and the NW clump at 
offset 8", 8" and 43.5 to 48.0 kms -1 may make another 
outflow), G8. 68—0. 37 is dominated by the simple, bipolar 
morphology discussed above which is wide-angled, espe- 
cially the red lobe. As shown in Figure [2] the CO outflow 
is well aligned with the shocked emission seen in IRAC 
bands presumably from the interaction of the outflow 
w ith the molecular envelop e . 

iCvganowski et all ([2009D find methanol masers are de- 
tected towards a large fraction of massive star form- 
ing regions associated with extended 4.5 p,m emission, so 
the association of the two in G8. 68— 0.37 is not surpris- 
ing. However, the resolution of the SMA observations 
is not sufficient to add to the ongoing debate whether 
these masers are asso ciated with outflows or discs (see 
IDe Buizer et al.ll2009l and references therein). 

4.3. Data Modelling 

4.3.1. Temperature and density profile modelling of 
continuum data 

Fol lowing lZhang et al.l ([20091 ) (see alsolTakahash i et al.l 
2009), we used the measured visibility amplitude of the 
SMA continuum data as a function of ui>-distance to de- 
termine the spatial density profile of G8. 68— 0.37. As- 
suming the density, p, and dust temperature, Tdust, scale 
as a power-law with radius (i.e. T^ustM oc r~ KT and 
p(r) oc r~ Kp ) and the dust is centrally heated, then the 
flux density from dust emission, F, is given by F oc J p 
TdustdS 1 , where S is the length along the line of sight. 
For k p + kt > 1, in the image domain the flux density 
scales as F oc r -( K p+ K r-i) j wn iJ e m the uv plane this 

corresponds to A oc S'™ P+KT , where A is the visibility 
amplitude and S uv is the mi-distance. Figure [5] shows 
the visibility amplitude as a function of wv-distance for 
the subcompact (triangles), compact (crosses) and ex- 
tended (circles) array configurations, scaled by the ob- 
serving frequency. The error bars show the uncertainty 
calculated from the variation in amplitude for the visi- 
bilities binned for each point. Assuming a temperature 
exponent kt = 0.33 (appropriate for the dust being cen- 
trally heated), the weighted least-squares fit to the data 
(shown as a line in Figure[6]) reveals k p = 1.8 ± 0.2. The 
intercept of the fit, corresponding to th e total flux , agree s 
well with the value of 3.4 Jy reported in lHill et all ([2005D . 

4.3.2. Radiative transfer modelling of NH3 data 

To determine a second, independent measurement of 
the spatial temperature and density profile we used 
radiat ive transfer modelling to fit the lLongmore et al.l 
(p007l) NH 3 (1,1), (4,4) and (5,5) data. The radiative 
transfer code, MOLLIE0, can deal with arbitrary 3D ge- 
ometries, but based on (i) the morphology of the single- 
dish sub-mm continuum emission, and (ii) analysis of the 
density profile in M4.3. It a spherically symmetric model 
was chosen as a good approximation, at least down to 

4 See IKetd (fT990T) for a description of the MOLLIE code, 
Keto et al.l (120041) for a description o f the line-fitting and 
Carolan et a l. (2009) and Kcto & Zhang (2010) for recent exam- 
ples of work using the code 
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spatial scales of ^0.05 pc probed by the extended con- 
figuration SMA observations. A model radius of 0.58 pc 
was chosen based on the extent of the single-dish sub- 
mm continuum emission. The temperature, Tr-, and 
density, p, were again parameterized with radius, r, as 
power laws: p(r) — pi/ 2 r~ Kp and Tj<-(r) = T X j 2 r~ KT , 
where pi/ 2 and T X j 2 are the density and temperature at 
the half-radius, 0.29 pc. An additional constant velocity 
component, AVnt, was included to model non-thermal 
support (e.g. from turbulence), resulting in a five pa- 
rameter model for the cloud. The NH3 to H2 abundance 
was fixed at 5 x 10~ 9 (derived from the observations of 
G8.68-0.37 in lPillai et aHl2007t ). This is lower than typi- 
cal abundances found in young massive star formation re- 
gions but increasing the abundance to several 10 -8 does 
not affect the general results of the paper. 

Models were constructed over a range of values for the 
five parameters (see Table [3|). Radiative transfer mod- 
elling was then used to generate synthetic data cubes 
for the NH 3 (1,1), (4,4) and (5,5) emission. These were 
then convolved with 2D Gaussian profiles at a spatial 
scale c orresponding to t h e reso lution (synthesized beam) 
of the iLongmore et al.l (|2007| ) observations. The syn- 
thetic spectra at each transition and position were fit to 
the observed spectra at three positions and reduced-^ 2 
values returned for the goodness-of-fit. The first posi- 
tion was located at the peak of the sub-mm continuum 
emission (aj 2 ooo =18:06:23.47, S J200 q = -21:37:7.6), 
then in radial steps of ~10" (the synthesized beam 
size) at aj 2 ooo =18:06:23.83, <5 ,72000 = —21:37:15.6 and 
aj20oo =18:06:24.33, S J20 oo = -21:37:21.6. While the 
NH3 integrated intensity emission shows a similar ex- 
tent to the single-dish continuum emission, the NH3 
morphology is extended along the outflow axis. As 
NH3 is known to be affected by outflow interactions 
(Zha ng et al.ll2007bl ). it is possible the NH 3 emission may 
be affected along this axis. Therefore, the radial direc- 
tion was chose n to be that perpendicular to the molecular 
outflow (see § 14.21 and Figures [2] and [5]) to minimize any 
potential contamination. 

Simulated annealing with 10,000 models was used to 
search through the 5D parameter space to minimize x 2 
and find the best- fit model. Figure H] shows the reduced- 
\ 2 values for the range of parameter space covered by one 
such run through the simulated annealing process. This 
method is inherently robust against becoming trapped 
in local, rather than global minima in parameter space. 
However, to determine the robustness of the best-fit 
model we ran the fitting 20 times with widely separated 
initial start values and increments. The results are simi- 
lar to those in Figure [8] Table [3] lists the range of param- 
eter space covered in the fitting proces s and the result- 
ing be st-fit values. Figure [7] shows the ILongmore et al.1 
(2007) NH 3 (1,1), (4,4) and (5,5) spectra from each of 
the three positions overlayed with the synthetic spectra 
at the same position from the best-fit model. 

4.3.3. CH 3 CN modelling 

Emission from the CH3CN J=12— 11 transition 
[~220.7GHz] is detected in the K = to 7 components^ 
in the 230 GHz compact-configuration observations. The 

5 Upper energy levels for K = to 7 components are 69, 76, 97, 
133, 183, 246, 325, 417 K, respectively 



emission is unresolved in the 3"x2" beam and at the 
peak of the mm continuum emission, encompassing the 
location of MM1 to MM3. Physical parameters of the 
gas at this spatial scale were estimated by x 2 mini- 
mization fitting of the spectrum, solving for the opti- 
cal depth assuming LTE conditio ns. With a collisional 
rate of 10" 8 cm 3 s _1 (|Greenl 119861 ) and Einstein A coef- 
ficient of 10 s _1 the transition has a nominal critical 
density of 10 4 cm~ 3 so the excitation should be collision- 
ally dominated in densities typical of massive star form- 
ing cores. The observed spectrum and resultant fit are 
shown in Figure [9l While the single-component model 
clearly fits the data well, the increasing linewidth and 
slight excess of emission in the high-K vs low-K compo- 
nents means the emission probably has more than one 
component, and the gas potentially has a temperature 
and density gradient. However, the signal-to-noise of 
the higher K-components is not sufficient to constrain 
more detailed modelling. The best-fit temperature in- 
cluding all K-components is 200 K which was robust to 
the 10% level under variation of the initial fitting es- 
timates and number of free parameters. Only fitting 
the K=0 to 4 components reduces the temperature to 
100 K but gives a much poorer result for the higher 
K components. The best-fit CH3CN co lumn density 
is 10 16 cm -2 . Gi ven typical abundances dWilner et al. 
[T99llZhang et al.lfl99cl (Hatchell et al.lll998t|Chen et al. 
120061: iZhang et al.ll2007at iLeurini et al.l 120071 ). the total 
H2 column density at this size scale is at least 10 24 cm -2 
and possibly significantly larger. For the above models, 
the filling factor implies a size of 0.8" for the CH3CN 
emission region. 

4.3.4. Model Consistency 

We now compare the results of the modelling out- 
lined in § 14.3.11 14.3.21 and 14.3.31 to check for consis- 
tency. At larger spatial scales, the results of the con- 
tinuum and NH3 modelling agree well, suggesting the 
simple spherical model with power-law temperature and 
density gradient is a good approximation at this spatial 
scale. However, the fit to the SMA continuum data in 
Figure [6] shows the assumption of spherical symmetry 
breaks down at angular scales below a few arcseconds - 
corresponding to the size at which the extended configu- 
ration resolves the emission into the three cores, MM1 to 
MM3. In the subsequent discussion section we refer to 
this radius of ^0.05 pc as the 'minimum reliable radius' 
for the power-law density and temperature approxima- 
tion. 

Given the poorer resolution (>8"), the observed NH3 
emission will be dominated by gas at much larger spa- 
tial scales. However, extrapolating the temperature and 
density gradient derived from the NH3 data down to 
0.8" (the best-fit size of the CH3CN emission) predicts a 
temperature of 100 K and density of 2xl0 7 cm~ 3 , which 
matches well with properties derived from the CH3CN 
data when only fitting the lower K-components. The 
excess emission and linewidth in the higher CH3CN K- 
components points to the density and temperature gra- 
dient continuing at even smaller spatial scales. In sum- 
mary, modelling of the three different data sets produces 
a coherent picture of the density and temperature profile 
of G8.68-0.37. 
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5. MASSES OF MM1 TO MM3 

We calculate the mass of MM1 to MM3 using the 
temperatures derived from the NH3 and CH3CN data 
- ^100 K at the radius encompassing all three com- 
ponents and likely increasing to ^200 K at smaller 
radii. Without independent temperature measurements 
of each component, the plausible mass range of MM1 
to MM3, assuming standard dust properties and follow- 
ing [Kau^annital] (|2008[) . is then 7-14M Q , 5—9 Mq 
and 2— 5M Q , respectively. This is a small fraction of 
the total mas s of 1500 M^ d erived from the single-dish 
observations (jHill et al.ll2010h . 

6. DISCUSSION 

We now return to the original focus of the paper to 
investigate the case for thermal feedback affecting frag- 
mentation in this young massive protocluster. 

We consider a simplistic scenario where, at t = io, 
the region initially had conditions similar to those of 
IRDCs: T - 10 - 2QK and n ~ 10 4 cuT 3 . A first 
generation of stars formed from this gas and began to 
heat up the surrounding environment. Some time later, 
at t = tfrag, the more massive cores we observe with 
the SMA formed from thermal fragmentation of the gas. 
The gas temperature and density at tf rag was Tf rag & 
nfrag- Assuming purely thermal fragmentation, this gas 
temperature and density set the mass (Mf rag ) and sep- 
aration (Af ra g) of the resulting fragments. Some further 
time later, at t — t now , we observed the region and 
have determined the present-day temperature [T now (r)] 
and density [/9 n ow( r )] as a function of radius from r = 0.05 
to 0.58 pc and the current mass (M now ) and separation 
(Anow) of cores MM1 to MM3 resolved on smaller scales. 

In the following sections we use these observational re- 
sults to infer the possible conditions at tf rag (i.e. Tf rag , 
/Oftag, Af rag and Mf ra J3). Given these inferred conditions, 
we then investigate the plausibility of different sources 
of feedback raising the gas temperature to the required 
Tfrag- In this way we aim to assess whether thermal 
feedback may have affected fragmentation. 

In § !6.1l we start with the assumption that the present- 
day density profile and core mass/separations are similar 
to those at the time of fragmentation - s o p now — pf rag , 
Anow = Af rag and M now =Mfra g . In § 16.21 we relax these 
assumptions. 

6.1. Thermal fragmentation assuming p now = /9f ra g, 
Anow = Afrag and M now = M[ rag 
6.1.1. Inferring Tfrag 

Assuming MM1 to MM3 have not moved since frag- 
mentation, their relative positions in Figure[5]then corre- 
spond to their fragmentation separations, Af rag , or Jeans 
length, projected on to a 2D plane. The measurement 
uncertainty in the projected separations is small - of or- 
der a few percent for the significance of the MM1 to 
MM3 detections (see e.g. lFomalontiri999l Eq 14-5), or 
calculated from the r.m.s phase noise, A(f>, and synthe- 
sized beam, 9b (positional uncertainty ~ A0 9b/2it). 
The major systematic uncertainty is the distance to the 

6 ^frag an d Mfrag & re equivalent to the Jeans length and Jeans 
mass, respectively. 



region, which we estimate to be ~20%. As the mea- 
sured separations are in projection, they underestimate 
the real separations. We try to quantify this underes- 
timation through geometrical arguments. Making the 
reasonable assumption that MM1 to MM3 lie within the 
volume encompassed by r < 0.05 pc (as opposed to lying 
at much larger radii and only seen in projection towards 
the center) , the largest underestimate will occur where a 
given pair lie at opposite sides (ie front vs back) of the 
r < 0.05 pc volume. Assuming MM1 lies near the center, 
the maximum separation is then 0.05 pc. The projected 
separations then underestimate the real separations by 
<40%. Any underestimate in separation will be reflected 
in an underestimate of Tf rag . 

Completeness is another potential observational bias 
in determining Af rag . If we are missing lower- mass frag- 
ments which fall below the detection limit, for example, 
the distance from MM1 to MM3 to the nearest fragment 
may be lower than the measured projected separations. 
However, from the sensitivity of the SMA observations 
(crms ~ 1 — 2mJy), any isolated fragments above ~1 Mq 
should be detected by their dust continuum emission. 
This 1 M Q limit is equal to the global Jean's mass for 
the physical conditions typical of cold IRDCs prior to 
the formation of any stars. A more in-depth discussion 
of the potential existence of a population of pre-existing 
but undetected lower mass stars is given in § 16.1.21 

A final concern in determining Af rag might be that 
MM1 to MM3 arc comprised of multiple, unresolved dust 
continuum peaks. For example, ^0.2" angular resolution 
observations towards similar regions at a similar distance 
and evolutionary stage resolve multiple (proto)stars with 
linea r separations down to 1700AU (|Longmore et al.l 
2006). However, such systems would correspond to mul- 
tiple stars forming within a single fragment, rather than 
reflecting the spatial distribution of fragments at tf rag - 

The uncertainties in the mass of MM1 to MM3, and 
hence Mf rag , are far larger than those of Af rag . In ad- 
dition to the large uncertainty introduced by the lack 
of independent temperature measurements for MM1 to 
MM3 (see §0, there are major systematic uncertainties 
in deriving gas masses from thermal dust emission - e.g. 
dust properties, distancf], gas-to-dust ratio etc. 

The density distribution is constrained at the larger 
spatial scales by the NH3 modelling (§ I4.3.2|) and at arc- 
second scales by the CH 3 CN modelling (§ I4.3.3|) . As 
discussed in § 14.3.41 these independently predict a den- 
sity of ~2xl0 7 cm~ 3 at the size scales of MM1 to MM3. 
This provides an upper limit to pf ra g as the rest of the 
volume containing MM1 to MM3 is likely to be at a lower 
density. As a more realistic estimate of pf ra g, the den- 
sity profile derived in § l4.3l predicts an average density of 
8xl0 6 cm~ 3 within a volume of r < 0.05 pc (the radius 
encompassing MM1 to MM3). 

Figure [TU] encapsulates the above uncertainties and il- 
lustrates how these affect the inferred fragmentation tem- 
perature. The solid lines show the expected Tf rag as a 
function of density for Jeans lengths of 0.05 and 0.03 pc, 
corresponding to upper and lower limit estimates of the 
MM1 to MM3 separations. The dotted lines show the 

7 The uncertainty in mass is proportional the uncertainty in 
distance, AD, squared, rather than linearly proportional to AD as 
is the case for the uncertainty in the measured separations. 
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expected Tf rag as a function of density for Jeans masses 
of 3 and 10 M Q , corresponding to reasonable estimates 
of the MM1 to MM3 masses. This shows that for the 
predicted densities, the fragmentation temperature must 
be at least 100 K. 

We now consider the plausibility of different sources 
of feedback raising the gas temperature to this required 

Tfrag ■ 

6.1.2. The case for heating by pre-existing, unobserved, 
embedded lower-mass stars 

From the sensitivity of the SMA observations (ctrms ~ 
1 — 2 mjy), any isolated dust continuum fragments above 
~1 M Q should be detected, but the SMA would not de- 
tect unembedded stars. The Spitzer observations (e.g. 
Figure [5]) would detect stars, but not if they were very 
deeply embedded. Due to spatial filtering, interferom- 
eters are only sensitive to density contrasts, so it is 
possible the observations are missing a population of 
uniformly-distributed, embedded lower mass stars. 

To assess the feasibility that pre-existing and unob- 
served stars are responsible for heating G8. 68— 0.37, we 
investigate what stellar population would be required to 
generate the required 100 K thermal fragmentation tem- 
perature. We take the simple approach of calculating the 
effect of heating by individual stars of different mass and 
luminosity to find the radius out to which they will heat 
the gas to 100 K. Given this radius, we then determine 
the minimum number of these stars uniformly spaced 
within the total spherical volume of radius 0.05 pc nec- 
essary to raise the temperature to 100 K. 

Assuming the stars and surrounding dust are in radia- 
tive equilibrium^ and the opacity can be approximated 
by a power-law of wavelength proportional to A~ p , a star 
with temperature, T*, and radius, i?*, will raise the dust 
to temperature, T d at a radius, r d (r<j >> -R*), given by 
r d ~ {Rj2)(T d /T„)-^+Py 2 (|Lamers fe Cassinellil Il99l 
Eq. 7.37). We adopt p = 1.5 based on the o pacities 
determined from the MRN (jMathis et al.l 119771 ) grain- 
size distribution between wavele ngths of 0.1/xm and 1mm 
(|Kruegel fc Siebenmorgen|[l994l ). To investigate the case 
of maximal heating, the central stars are assumed to 
have reached the main sequence. This provides an up- 
per limit to the stellar effective temperature because 
protostars on their Hayashi tracks are generally cooler 
(|Coxil2000L Table 15.14). Figure [TT] shows the resulting 
dust temperature as a function of radius from stars of 
mass = 0.8, 1.2, 1.8, 2.5, 5 and 10 M©, overlayed with the 
best-fit model profile determined in § 14.31 Main-sequence 
stars of 0.8, 1.2 and 1.8 M Q will heat dust to >100K at 
radii < 3.5xl0~ 4 , 1.0xl0~ 3 and 3.5xl0" 3 pc, respec- 
tively. To heat the volume of radius 0.05 pc to 100 K in 
this way therefore requires uniformly distributing these 
stars by the above separations throughout the volume. 
The minimum number of stars this implies (2.9 xlO 6 , 
1.25xl0 5 and 2.9xl0 3 ) within a radius of 0.05 pc leads 
to total stellar masses far in excess of the total mass of 
gas available. We therefore conclude the number of stars 
required to heat the gas is unfeasibly large. 

At the earliest evolutionary stages, the luminosity from 
low-mass stars can be dominated by accretion rather 

8 i.e. each radius is considered as a geometrically thin and 
optically-thick shell with no opacity between it and the star. 



than nuclear burning. We now consider the affect of heat- 
ing due to accretion luminosity. As most of the infalling 
particles' potential energy will be released at radii close 
to the protostar, the luminosity can be approximated as 
arising from a point-source. Following similar assump- 
tions about the dust properties etc. as before, the dust 
temperatu re as a function of radius can be solved analyt- 
ically (e.g. iScoville fc Kwanlll976l: iGarav fe Lizanolll999l 
Eq 11). The left-hand panel of Figure [T2l shows the dust 
temperature as a function of radius for central heating 
sources with luminosities from 1 — 10 5 L Q . This range was 
chosen to include the expected luminosity from accret- 
ing low-mass stars (which are dominated by accretion- 
luminosity) through to O-stars dominated by stellar lu- 
minosity. The right-hand panel of Figure [12] shows the 
radius at which central heating sources of a given lu- 
minosity will heat the surrounding gas up to >100K. 
Following the same line of argument as before, we can 
investigate how many sources are required to heat the 
100 M volume of radius 0.05 pc (the mass determined 
from the modelling that lies within the radius encompass- 
ing MM1 to MM3) to 100 K as observed in G8.68-0.37. 
iFroebrichl (|2005l ) show that isolated Class cores (which 
as the youngest protostars should have the highest accre- 
tion rates) typically have luminosities <1OL0. Again, an 
unfeasible number of low mass stars are needed to heat 
the gas to the required temperature, even allowing for the 
additional heating from accretion luminosity. Assuming 
the higher infall rates in a massive protocluster environ- 
ment give rise to larger accretion luminosities does not 
solve the problem either. The total luminosity required 
would be larger than the present-day measured bolomet- 
ric luminosity of 1 .9xl0 4 LfD. 

Comparing the lLamers fc Cassinclli (1999) tempera- 
ture profiles to those of IGarav fc Lizanol (17999 Fl we 
note the former are factors of a few higher. The 
IGarav fc Lizanol (I1999T) formulation app ears closer to ob- 
served profiles (e.g. lKeto fc Z hang 2010, and the gas tem- 
peratures derive d from our own obse r vation s (§ 14. 3|1 ) so 
we consider the lLamers fc Cassinellil (|1999l ) formulation 
(e.g. Figure ITT1) as an upper limit to the temperature. 

The discrepancy between these two methods does not 
affect the conclusion from the above analysis. A large 
population of low mass stars would not be able to heat 
the gas to the 100 K thermal fragmentation tempera- 
ture required to reproduce the observed separations and 
masses of cor es MM1 to MM3. This is consistent with 
the models of lOffner et al.l (|2009l ) which show that low- 
mass stars have a negligible affect on the overall cloud 
heating. Repeating the above calculations but imposing 
more realistic limits - that the total stellar mass is not 
larger than the available gas mass and the total luminos- 
ity is not larger than the measured bolometric luminosity 
- gives an upper limit to the temperature provided by 
low-mass stars of ~40K. 

6.1.3. The case for heating by an early-B star 

The calculations in the previous section show that low- 
mass stars would not be able to heat the gas within a 
radius of 0.05 pc to the temperature of 100 K required 

9 We compared the two analytic solutions by taking the known 
radii and effective temperatures of ZAMS stars of a given luminos- 
ity 
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for the observed cores (MM1, MM2 and MM3) to have 
formed by thermal fragmentation. The solid line rep- 
resenting the best-fit model temperature profile in Fig- 
ure QT] shows that this heating could be attributed to a 
single early B-type star rather than a larger number of 
low-mass stars. However, invoking an early B-type star 
as a source of thermal feedback to raise the gas tem- 
perature and suppress fragmentation (allowing cores as 
massive as MM1, MM2 and MM3 to form) is unsatisfac- 
tory as it invokes a circular argument - it does not solve 
the problem of what initially raised the temperature to 
suppress fragmentation into lower mass fragments which 
would have allowed the B star to form in the first place. 

6.2. Thermal fragmentation relaxing the assumptions 
that p now = Pfrag, A n0 w = Af rag and M now = Mf rag 

Having concluded it is not possible for thermal feed- 
back to have raised the gas temperature to >100K if 
Pnow = Pfrag, A now = A frag and M now = Aff rag , we con- 
sider how changes in the gas conditions from if rag to the 
present-day could have affected the resulting fragments 
and look for observational predictions that these changes 
would imply. 

As shown in Figure I1Q[ the required fragmentation 
temperature would be lower than 100 K if the density 
at the time of fragmentation were lower than the present 
day - ie if p n0 w > Pfrag • I n § 16.2.11 we investigate the ob- 
servational evidence for the large-scale infall of gas that 
this would imply. In § !6.2.2l we consider how a lower den- 
sity of gas at £f rag would have affected the separations of 
MM1 to MM3 over time. 

6.2.1. The case for large-scale infall 

As discussed i n § 1 1.11 the molecular line profiles 
in IPurcell et al.l (|2009T > provide observational evidence 
for large-scale infall of gas in G8.6 8 — 0.37 . Follow- 
ing the calculations in Wals h et al.l ((2006), the in- 
fall rate for G8.68-0.37 inferred from the HCO+ and 
HCN spectra is ~10~ 4 M yr" 1 . This infall rate is 
significantly larger than typically measured towards 
lower mass star-forming cores but comparable to ob- 
servatio ns of other young high mass star formation 
regions flWu fc Evansi [ 2003; Pe retto et all 120061 120071: 
Wals h~etar.ll2006t iKeto fc KlaasserJ I2008t Ide Wit et al.l 
20091) . 

Are these infall motions suggested by the molecular 
line observations consistent with the velocities expected 
in gravitational collapse? We estimate the dynamical 
state of the region by comparing the contribution from 
thermal and non-thermal support as a function of radius 

to th e virial velocity, Kiriai = (tt) 1/2 (|Stahler fc Pallal 
2005, Eq 3.20) - a measure of how much kinetic energy 
is required to balance the gravitational potential. To 
convert from the measured one-dimensional velocity to 
the three-dimens i onal r oot mean square velocity we use 
iRohlfs fc Wilson! (|2004f ). Eq 12.72. The virial velocity 
in Figure [T3] is then calculated at each radius using the 
enclosed mass (derived from the model parameters in 
§ 14. 3|) inside that rad ius. The range o f uncertainty in 
virial velocity (see e.g. rElmcgrcer] |1989f) is illustrated by 
the hatched area between the dash-dot-dot lines. Fig- 
ure [13] shows, i) as expected the linewidth is dominated 
by the non-thermal contribution, and ii) even the higher 



temperatures towards the center are insufficient to reach 
the required kinetic energy support for the cloud to be in 
equilibrium. Although undoubtedly oversimplifying the 
gas dynamical state (both support from magnetic fields 
and the surface terms in the virial equation are ignored, 
for example), we find that the inward motions are gen- 
erally consistent with gravitational contraction. 

6.2.2. Could the cores have been more widely separated in 

the past? 

If the entire star forming cloud is globally contract- 
ing, then the density of the gas from which MM1 to 
MM3 fragmented could have been lower than measured 
today. The cores could initially have formed at larger 
separations and moved closer over time to reach their 
current locations. To test the feasibility of this scenario 
we can estimate how far these cores may have moved 
over a given time period and predict their resulting rel- 
ative velocities. Assuming the fragmentation conditions 
were similar to those of IRDCs (10 5 cm~ 3 & 10 - 20 K), 
which are thought to be massive star formation regions at 
the earliest evolutionary stages, the Jeans length would 
have been ^0.07 pc. To reach their current separations 
of ~0.048pc, the cores would have had to move ~0.01 pc 
towards each other. 

We use the free-fall time to derive a lower limit on how 
long the cores may have been moving since they formed. 

rrm ( <t \ 1/2 

Figure [Ml shows the free-fall time, tff = ( 32 q p J > 

of the gas in G8. 68— 0.37 as a function of radius, where 
p av is the average density within each volume calculated 
from the enclosed mass within that radius. Taking 10 yr 
(the free-fall time at 0.05 pc) as a lower limit of the proto- 
cluster age and following I Walsh et al.l (|2004h . cores mov- 
ing 0.01 pc would reach a velocity of 2.6 kms -1 . Future 
molecular line observations with sufficient angular reso- 
lution to identify MM1 to MM3 should easily be able to 
discern such large relative velocities. The inferred proper 
motions of ~60/x"/yr are too small to detect. 

6.3. Considering non-thermal fragmentation 

In the preceding sections we have used Jeans analy- 
sis to determine the mass and separation of core frag- 
ments, given the temperature and density of the sur- 
rounding environment. While the assumption of purely 
thermal support is undoubtedly an oversimplification - 
both turbulence and magnetic fields are also likely impor- 
tant as support mechani sms at the earliest stages of mas- 
sive star formation (e.g. [Zhang et al"1l2009t iGirart et al.l 
l2009t iTang et al.l I2009al lb n2010l) - the relative impor- 
tance of the different support mechanisms is not well un- 
derstood at these spatial scales and for objects at such an 
early evolutionary stage. Without a direct measurement 
of the magnetic field towards G8. 68— 0.37 we make the 
assumption of equipartition between thermal, turbulent 
and magnetic-field support and assess how this affects the 
fragmentation temperature determined in § 16.1.11 Re- 
turning to Fig 110) we used a lower-limit fragmentation 
temperature of 100 K in the analysis of § 16.11 If tem- 
perature, magnetic fields and turbulence are in equipar- 
tition the total energy is equivalent to a temperature of 
300 K. Fig [T0l shows that 300 K is in fact a more realis- 
tic fragmentation temperature than the 100 K previously 
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assumed. We therefore conclude that the results derived 
assuming purely thermal fragmentation are robust when 
including support from turbulence and magnetic helds. 

Finally, we consider the assum ption in the Jeans anal - 
ysis that the gas is isothermal. iKrumholz et al.1 ()2007f ) 
find in their numerical models that at high column den- 
sities, where the gas becomes optically-thick, the equa- 
tion of state deviates from isothermal. The result is that 
fragmentation is suppressed at a lower temperature than 
would be expected from Jeans analysis. While this is 
difficult to test with observations of a single region, in 
future work we will look to compare the mass distribu- 
tion of cores above and below the critical column density 
of lgcm~ 2 , at which this affect is pr edicted to become 
important (jKrumholz fe McKeel [20081 ). 

6.4. Will the G8.68- 0.37 protocluster form O stars? 

In order to form an O star through direct collapse, a 
fragment would need to have a mass of at least the stel- 
lar mass times by the star formation efficiency - ie of 
order 100 M q of gas, given typi cal star formation effi- 
ciencies (e.g. lLada fc La da 2003). No matter what the 
fragmentation history of the region, none of the observed 
MM1 to MM3 cores has this much mass. Any feedback 
(if it exists/existed) was therefore not sufficient to create 
a fragment that can collapse to form an O star. In or- 
der to form an O star via accretion, these existing cores 
need to obtain additional mass from elsewhere. We now 
investigate whether G8. 68— 0.37 is likely to form O stars 
and, if so, how this might occur. 

The total mass within r < 0.05 pc (encompassing MM1 
to MM3) predicted from the best-fit model density profile 
is ~100M Q . MM1 to MM3 account for roughly a third of 
the mass, leaving a large reservoir of material inside this 
relatively small volume. If the measured infall rate (see 
§ 16.2. lj) can be sustained for a few 10 5 yr and large-scale 
infall can feed mass to smal ler spatial scales (as observe d 
towards G20.08-0.14N bv IGalvan-Madrid et al.l [2009h . 
it seems plausible that MM1 to MM3 may end up as 
O stars via continued accretion. However, without in- 
formation of the gas velocities at r < 0.05 pc, it is not 
clear how MM1 to MM3 are coupled to the gas at this 
size scale and whether or not they will continue to gain 
mass. Several theoretical scenarios are postulated: con- 
tinued accreti on through a common rotating envelop e 
or 'disk' (e.g. iKetol [2003L [20071 IKrumholz et all [20(11 : 
acc retion from initial ly unbound cluster-scale gas (see 
e.g. lSmith et al.ll2009| ). or fragmentation-induced starva- 
tion via grayitatiq nal instabilities in the accretion flow 
(|Peters et al.ll2010fl . 

Considering the first scenario - that MM1 to MM3 
may be fed by accretion through a common rotating en- 
velope or disk - we do not detect molecular-line emission 
in the highest angular resolution observations so can not 
search for kinematic signatures of rotation at this spa- 
tial scale. However, it is interesting to note that MM1 to 
MM3 are oriented roughly in a plane per pendicular to the 
outflo w. This is similar to the results of iFallschecr et al.l 
(2009) towards 18223-3 - another young high mass pro- 
tocluster, with similar global mass, infall and outflow 
rates. They find the mm continuum emission peaks are 
aligned perpendicular to the outflow and encompassed 
by a flattened, rotating entity of inward spiralling molec- 
ular gas. Although emission from a rotating envelope 



of 28,000 AU would have been resolved in the observa- 
tions of G8. 68— 0.37, it is not seen. But such a large 
rotating envelope is not necessarily expected. A recent 
large survey to find high mass accretion disk candidates 
(jBeuther et al.ll2~009h found diameters of the rotating en- 
tities to be ~10,000AU, smaller than the unresolved 
bright molecular line emission in the SMA observations 
in this work. 

The second scenario is similar to that described in 
§ 16.2.21 where the cores initially form with smaller mass 
and at larger radii. Over time they move closer under 
the influence of the global gravitational potential. In this 
picture, cores grow by accreting the diffuse surrounding 
unbound cluster-scale gas. With this in mind, it is in- 
teresting that the most massive observed mm continuum 
core, MM1, is found at the center of the gravitational po- 
tential (as predicted) where the infall would be greatest. 
However, other than searching for the relative velocity of 
the cores (see § I6.2.2[) it is difficult to test these predic- 
tions with observations of a single region. We are cur- 
rently in the process of extending this analysis to larger 
samples of massive star formation regions (e.g. from the 
NH3 cores detecte d in the HOPS Galactic Plane survey 
IWalsh et al.ll2008h . 

In summary, we conclude that G8. 68— 0.37 may still 
proceed to form O stars, with large scale infall and con- 
tinued accretion at smaller scales feeding the observed 
mm cores. Observations to resolve the gas kinematics 
at r < 0.05 pc are required to determine i) how the mm 
cores and remaining gas are coupled at this size scale 
and ii) the mechanism through which accretion is taking 
place. 

7. CONCLUSIONS 

Combining multiple-configuration SMA continuum ob- 
servations with modelling of the NH3 and CH3CN 
emission, we have determined the present-day i) tem- 
perature/density as a function of radius and, ii) 
mass/separation of cores, in the massive protocluster 
G8. 68— 0.37. We then used this investigate whether feed- 
back from low mass stars can raise the protocluster gas 
temperature sufficiently to delay thermal fragmentation, 
allowing massive stars to form through direct collapse of 
high-mass fragments. 

From radii of 0.58 pc (1.2xl0 5 AU) down to 0.05 pc 
(10 4 AU) we find the data are well fit with the region 
having a power law temperature and density of the form 
Toe r~ - 35 and p oc r ~ 2m . At r < 0.05 pc the assump- 
tion of spherical symmetry breaks down and the tem- 
perature/density range are calculated to be 100— 200K 
and 10 6 — 10 7 cm -3 , respectively. Within this radius of 
0.05 pc, the SMA resolves the 1.2 mm continuum emis- 
sion into 3 cores, MM1 to MM3, with separations of 
6200 AU and 9700 AU and masses 10±3 M Q , 7±2 M Q and 
4±2M Q , respectively. 

From the region's observed properties, we infer the 
conditions at the time the cores formed by fragmenta- 
tion. Assuming the measured separations of MM1 to 
MM3 and the average density encompassing these cores 
(r < 0.05 pc) are representative of the physical conditions 
at fragmentation, this implies a thermal fragmentation 
temperature of at least 100 K. We rule out a population 
of low mass stars being able to provide this heating - an 
unfeasibly large number are required and the measured 
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bolometric luminosity is orders of magnitude too small 
even allowing for heating due to trapped accretion lu- 
minosity. The heating could instead be provided by the 
equivalent of a 5 M© ZAMS star but this is an unsatis- 
factory source of thermal feedback - it invokes a circular 
argument requiring the star to have raised the tempera- 
ture to suppress fragmentation, before it formed. Alter- 
natively, the required fragmentation temperature could 
be lower if the region were initially at a lower density, for 
example if the cores formed farther apart or the region 
was undergoing global infall. 

Whatever the fragmentation history, none of the ob- 
served cores has sufficient raw material to form an O 
star through direct collapse. Even if feedback may have 
suppressed fragmentation, it was not sufficient to halt it 
to this extent. If G8. 68— 0.37 is destined to form O stars, 
the observed cores must obtain additional mass from out- 
side their observationally derived boundaries. The ob- 
servations suggest that the cores in this protocluster are 
being fed via global infall from the very massive reser- 



voir (^1500 M Q ) of gas within which the protocluster is 
embedded. 
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Overview of the observing setups. The first two columns show the array configuration and observing frequencies. The 
sub-compact, compact and extended array configurations have baselines ranging from 31 - 89 m, 61 - 192 m and 82 - 247 m, 
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Properties of the mm continuum sources detected with the SMA in the extended array configuration. 
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TABLE 3 

NH3 RADIATIVE TRANSFER MODELLING PARAMETERS. A SPHERICALLY-SYMMETRIC MODEL OF RADIUS 0.58 PC WAS CHOSEN BASED ON THE 
EXTENT/MORPHOLOGY OF THE NH3(1,1) EMISSION. THE TEMPERATURE, T K , AND DENSITY, p, WERE PARAMETERIZED WITH RADIUS, r, AS 

POWER LAWS: p(r) = p 1 l 2 r _re P AND T K (r) =Ti/2 r~ KT , WHERE p 1 l 2 AND T1/2 ARE THE DENSITY AND TEMPERATURE AT THE HALF 
RADIUS POINT, 0.29 PC. An ADDITIONAL CONSTANT VELOCITY COMPONENT, AVnT, WAS INCLUDED TO MODEL NON-THERMAL SUPPORT (E.G. 
FROM TURBULENCE), RESULTING IN A 5 PARAMETER MODEL FOR THE CLOUD. THE SECOND AND THIRD COLUMNS SHOW THE RANGE OF 
PARAMETER SPACE THE MODELS SEARCHED THROUGH. THE FINAL COLUMN GIVES THE BEST FIT MODEL TO THE L07A DATA AFTER 
MINIMISATION OF 10000 MODELS USING SIMULATED ANNEALING TO SEARCH TH ROUGH THE 5D PARAMETER SPACE. FULL DETAILS OF THE 

MODELLING ARE GIVEN IN § 14.3.21 



Parameter 


Search range 


Best-fit 




Min 


Max 


Value 


log (P1/2 [cm- 3 ]) 


3.0 


8.0 


4.98 




0.5 


4.0 


2.08 


T1/2 [K] 


10.0 


400.0 


41.5 




0.1 


2.5 


0.35 


AV NT [kms" 1 ] 


0.5 


3.0 


1.79 
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Fig. 1.— The IRAS 18032-2137 star forming complex . A GLIMPSE 3-colour image at 3. 6pm, 4.5pm and 8.0pm is overlayed with SCUBA 
850 pm continuum emission contours (Hillct al. 200(J. The complex compri ses 3 distinct regions: i) the stellar cluster, BDS2003 — 3 
IIBica et al l 12003]) . with the approximat e extent reported by IBica et al.l 12003) indicated by the circle; ii ) the ultra-compad 1111 region 
G8.67 — 0.36 (Wood & Churchwcll 1983, cm-continuum position illustrate d by box ); iii) a massive star forming core prior to formation of 
an HII region G8.68- 0.37. H 2 Q , CH 3 OH and OH maser emission l|Hofner fc ChurchwelllH996l : IWalsh et al.|[l998T : IForster k, Caswell|[T989l : 
ICaswell|[l998l:IVartts et al.ll2000l) are s hown as blue, black and red crosses, respectively. The line towards the lower left gives the physical 
scale for the source distance of 4.8 kpc HPurcell et al. l l2006T ). 
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18 h 06 m 24.5 s 24.0 s 23.5 s 23.0 s 22.5 s 



Right Ascension (J2000) 

Fig. 2. — 3-color GLIMPSE image at 3.6, 4.5 and 8.0 fim overlayed with (i) the 217 GHz SMA extended-configuration continuum image 
[grey (— 3<r) and black (3, 6, 9a, ...) contours: a = 2.6mjy] and (ii) the 12 CO(2— »1) emission integrated from 29.3-36.3 km/s [blue 
contours] and 45-74 km/s [red contours]. The blue and red lobes of the molecular outflow lie along a similar axis as the extended 4.5 lira 
emission - a well-known tracer of gas shocked in outflows. The 3 SMA 1 mm emission peaks are labelled MM1 to MM3 in order of peak 
1 mm continuum emission flux. The brightest component, MM1, is found at the center of the outflow lobes and appears the best candidate 
for the outflow driving source. The filled black cirle in the lower-left corner gives the synthesized beam of the SMA extended configuration 
observations. The bar in the upper-right corner illustrates the physical scale for the source distance of 4.8 kpc. 
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Fig. 3. — Lower sideband spectrum of the 230 GHz, compact-configuration SMA data taken at the position of peak continuum emission 
and smoothed to a resolution of 1.75 km/s. Transitions with 2 consecutive channels >3cr (cr = 92 mjy/channel) are labelled above the 
spectrum at frequencies (corrected for the source Vlsr of 37.2kms _1 ) taken from the NIST catalog (Lovas 2004). In addition to the 
13 CO (2— ¥l) and 12 CO (2— ¥l) emission, the spectra show several more complex molecules (e.g. CH3OH & CH3CN) confirming the source 
is warm enough for these molecules to have evaporated off the dust grains and into the gas phase. However, the spectra do not display the 
rich inventory seen towards hot molecular cores suggesting this region is an intermediate evolutionary stage between the cold and hot core 
stages. 
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Fig. 4. — Upper sideband spectrum of the 230 GHz, compact-configuration SMA data taken at the position of peak continuum emission 
and smoothed to a resolution of 1.75 km/s. Transitions with 2 consecutive channels >3cr (cr = 94mjy/channel) are labelled above the 
spectrum at frequencies (corrected for the source Vlsr of 37.2kms _1 ) taken from the NIST catalog (Lovas 2004). In addition to the 
13 CO (2— ¥l) and 12 CO (2— ¥l) emission, the spectra show several more complex molecules (e.g. CH3OH & CH3CN) confirming the source 
is warm enough for these molecules to have evaporated off the dust grains and into the gas phase. However, the spectra do not display the 
rich inventory seen towards hot molecular cores suggesting this region is an intermediate evolutionary stage between the cold and hot core 
stages. 
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Fig. 5. — 12 CO(2— >1) channel map towards G8.68 — 0.37. The Vlsr ( m kms — 1 ) of each frame is shown in the top-left of each panel 
and the synthesised beam size is shown as a filled ellipse in the lower left corner of each frame. The center of the coordinate system and 
location of the cross show the methanol maser position at aj2Q00 = 18:06:23.47, <5j2000 = —21:37:10.6. The contours start at ±5<r and 
increment in 5<r intervals (<r = 140 mjy/beam). For reference, the three stars show the location of MM1, MM2 and MM3. 




Fig. 6. — UV-distance vs visibility amplitude for the 1mm continuum emission observed with the SMA. Visibilities from the compact, 
subcompact and extended array configurations are shown as triangles, squares and circles respectively. The error bars show the uncertainty 
calculated from the variation in amplitude for the visibilities binned for each p oint. T he solid line shows the weighted least-squares fit to 
the data - assuming Toe r -0 - 33 as expected for an internally heated source (see § 14.3.11) the best-fit density profile is given by p oc r -l-8±0.2 
The dashed lines show the corresponding angular scale at 10" (left), 3" (center) and 1" (right). While the assumption of spherical symmetry 
appears a good fit to the data at large spatial scales, it is clear this assumption breaks down at < 3", where the extended configuration 
observations resolve MM1, MM2 and MM3. 
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Fig. 7. — Results of the NH3 radiative transfer modelling outlined in § 14.3,21 Longmorc ct al. (2007) spectra are shown in black, overlayed 
with the best-fit model spectra in red. Each column shows spectra taken from a single position on the sky. The top, middle and bottom 
rows are the NH3(1,1), (4,4) and (5,5) transitions, respectively. Th e posi tions w ere selected to be spaced by roughly one synthesised beam 
in the direction perpendicular to the molecular outflow (see § 14.21 and § I4,3.2[l to minimize potential contamination from gas affected by 
the outflow. 
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Fig. 8. — Reduced-x 2 values from the simulated annealing process to find the best-fit model to the Longmore ct al. (2007) data (see 
§ 14.3- 2^ . Each of the crosses shows the physical parameters (density, density exponent, temperature, temperature exponent and turbulent 
linewidth) from one of the 10,000 models, and the resulting reduced-x 2 value for that particular model. 
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220.75 220.70 220.65 220.60 220.55 220.50 
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[G. 9. — CH3CN (J=12— 11) spectrum from the 230 GHz, compact-confi guratio n SM A data at the position of peak continuum emission 
,ck] overlayed with the single-component LTE model fit as outlined in § 14.3.31 The K-components are labelled above the spectrum. 
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Fig. 10. — Encapsulating the observational uncertainties in the measured separations and mass of MM1 to MM3 to illustrate how these 
affect the inferred fragmentation temperature. The solid lines show the expected Tf rag as a function of density for Jeans lengths of 0.05 and 
0.03 pc, corresponding to upper and lower limit estimates of the MM1 to MM3 separations. The dotted lines show the expected Tf rag as a 
function of density for Jeans masses of 3 and 10 Mq, corresponding to reasonable estimates of the MM1 to MM3 masses. The dot-dashed 
lines illustrate the most extreme mass range given the observational uncertainties. The vertical dashed line shows the average density within 
a volume of r < 0.05 pc (the radius encompassing MM1 to MM3) of 7.8xl0 6 cm~ 3 as a reasonable estimate for the fragmentation density. 
The shaded region illustrating the observational limits shows the fragmentation temperature is at least 100 K (the horizontal dashed line). 
If the density were lower in the past (for example if the region were undergoing global infall) the fragmentation temperature would also be 
lower. 



Fragmentation & protostellar heating in a massive protocluster 



21 



o 




Core Radius 



10 3 0.01 0.1 1 

Radius (pc) 



Fig. 11. — Te mper ature as a function of distance from stars of various mass. The solid line shows the temperature of the best-fit model 
determined in § 14.31 fToc r~ 0,35 and p oc r -2 08 ). The right and central vertical dashed lines show the size-scale over which this profile 
is reliable - from the core radius down to 0.05 pc (where the mm-continuum emission is resolved into multiple components). The left 
vertical dashed line gives the highest resolution of the SMA observations. The dotted lines show the analytical relations hip between dust 
temperature and radius, assuming the dust is being heated and is in radiative equilibrium with a central star (Lamcrs & Cassinclli 1999, 
Eq. 7.37). To ensure an upper-limit to the stellar temperature, the stars are assumed to have reached the main-sequence (Cox 2000, Table 
15.14). Temperature profiles are shown for stars of mass =0.8,1.2,1.8,2.5,5 and 10 Mq. 
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Fig. 12. — [Left] Dust temperature as a function of radius for central heating sources of different luminosity. The right and left vertical 
dashed lines show the radius where the mm-continuum emission is resolved into multiple components, and the highest resolution of the 
SMA observations, respectively. The dotted lines show the analytical relationship between dust temperatur e and radius, assum ing the 
dust is being heated and is in radiative equilibri um with a central power ing source of the listed luminosity (Scovillc & Kwan 1976, Eq 9). 
[Right] Using the same analytical description in Scovillc & Kwan (1976), the solid line shows the source luminosity required to heat the 
gas to >100K at the distance given on the x-axis. 




Fig. 13. — Comparison of non-thermal and thermal support in G8. 68— 0.37. The horizontal da shed li ne, indicated as Vnt> shows the non- 
thermal contribution to the line- width determined from the NH3 radiative transfer modelling in S 14.3.21 The dot-dash line, labelled Vthermal, 
shows the expected the rmal contribution to the NH3 line-width as a function of radius, determined from the temperature distribution of 
the best-fit model in § 14.31 The solid line, Vihermal+NTi shows Vjjt an d ^thermal added together in quadrature giving the expected 

total linewidth. Finally, the dash-dot-dot lines and enclosed hatched area gives the expected range in virial velocity, V^irial = ( ) 

— a measure of how much kinetic energy is required to balance the gravitational potential of a cloud. The right and left vertical dot- 
dashed line shows the minimum reliable radius for the assumed power-law density distribution and the maximum resolution of the SMA 
observations, respectively. This figure shows, i) as expected the linewidth is dominated by the non-thermal contribution, and ii) even the 
higher temperatures towards the center are insufficient to reach the required kinetic energy support for the cloud to be in equilibrium. 
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32 ^ — J , as a function of radius where p av is the average density enclosed 

within a given radius determined from the enclosed mass within the radius on the horizontal axis. The right and left vertical dot-dashed line 
shows the minimum reliable radius for the assumed power-law d ensit y distribution and the maximum resolution of the SMA observations, 
respectively. The dynamical time of the outflow determined in ^ 14.21 is shown as a dashed horizontal line. The free-fall time at the radius 
at which the core is observed to fragment is similar to the dynamical time of the outflow. 



